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Abstract 

We consider the time evolution of systems in which a spatially homogeneous 
scalar field is coupled to fermions. The quantum back-reaction is taken into 
account in one-loop approximation. We set up the basic equations and their 
renormalization in a form suitable for numerical computations. The initial 
singularities appearing in the renormalized equations are removed by a Bo- 
goliubov transformation. The equations are then generalized to those in a 
spatially flat Friedmann- Robert son- Walker universe. We have implemented 
the Minkowski space equations numerically and present results for the time 
evolution with various parameter sets. We find that fermion fluctuations are 
not in general as ineffective as previously assumed, but show interesting fea- 
tures which should be studied further. In an especially interesting example 
we find that fermionic fluctuations can "catalyze" the evolution of bosonic 
fluctuations. 
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1 Introduction 



In recent years the study of nonequilibrium dynamics in quantum field theory 
has received much attention. Quantum fields out of equilibrium can play an 
essential role, e.g., in cosmology [1-16], in the quark-gluon plasma [17-25] or 
during phase transitions in solid state physics [26]. 

While the formalism of nonequilibrium dynamics in quantum field the- 
ory has been established long ago p7l El], real time simulations for realistic 



systems have been developed only recently. Numerical simulations of the evo- 
lution equations have been studied by various authors. The general features 
are similar: The quantum back-reaction cannot be described by Markovian 
friction terms. The relaxation of the classical field amplitude either shuts off 
or is powerlike, the quantum ensembles generated are characterized by para- 
metric resonance bands, the full development of the resonance being sup- 
pressed by the quantum back-reaction. The early- and late-time behavior 
has been analyzed analytically for the one-loop and large-N approximations 



29], |0|. 

Most of the numerical simulations have been performed with a scalar 
classical field coupled to scalar quantum fields. The quantum back-reaction 
of fermion fields on a classical scalar field has received little attention up to 



now. In j31| it has been stated, on the basis of some numerical evidence, 
that at least for large field amplitudes Pauli-blo eking would make fermions 
ineffective for dissipation and damping of the classical field. In another re- 
cent publication |32|] the leading orders in perturbation theory have been 
evaluated, so far without numerical computations. The interaction between 
a classical electric field and fermionic fluctuations has been considered in 



P"7J (see also |19[) as a model for qq production in the quark-gluon plasma. 
There the evolution of the system in quantum field theory was compared to 
the evolution using the Boltzmann-Vlasov equation. 

In this paper we reconsider fermionic fluctuations coupled to a scalar 
field, with respect to both formal aspects and numerical simulations. We use 



a formalism developed by us recently |33|, |3J] to formulate the renormalized 
equations of motion in a form which is suitable for numerical computation but 
satisfies at the same time the usual requirements of renormalized quantum 
field theory. The renormalization scheme is covariant and independent of the 
initial conditions. After renormalization we find the equations to be singular 



at t — 0, a phenomenon known as Stueckelberg singularities [35, 36|. As in the 



scalar case we have studied previously |3~7 1 , these singularities can be removed 



by a Bogoliubov transformation of the initial fermionic quantum state. We 
then generalize the equations to those in a flat FRW universe. We finally 
formulate the linearized equations of motion, in order to be able to compare 
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with the full quantum evolution. We have implemented numerically the 
formalism developed in this article. Numerical results for various parameter 
sets, and various aspects of these results will be presented and discussed in 
sections 8 and in the conclusions. 

The plan of this article is as follows: In section 2 we formulate the basic 
relations and the equation of motion for a scalar field coupled to fermions; in 
section 3 we present the energy momentum tensor and discuss the fermion 
number; the renormalization of the equation of motion and of the energy 
momentum tensor is developed in section 4; in section 5 we derive the Bo- 
goliubov transformation which removes the initial singularities; the extension 
to a conformally flat Friedmann- Robertson- Walker is derived in section 6; in 
section 7 we discuss the linearized equations of motion; the results of our 
numerical simulations are discussed in section 8; conclusions are given in 
section 9. 



2 Basic Relations and equation of motion 

We study a model consisting of a scalar " inflaton" field $ coupled to a 
spin 1/2 field ip by a Yukawa interaction. We do not introduce a genuine 
mass term for the fermion field, it acquires a time-dependent mass via the 
Yukawa coupling. We introduce a A$ 4 self-interaction, but do not consider 
the case of spontaneous symmetry breaking and do not consider the quan- 
tum fluctuations of the scalar field itself. This framework is sufficiently gen- 
eral for discussing renormalization and the typical effects introduced by the 
back-reaction of the fermion fields. It can be easily generalized to include 
a fermion mass, bosonic fluctuations, and spontaneous symmetry breaking. 
The Lagrangian density is given by 

£ = -d^d»<5> - -M 2 § 2 - + - g$)rp , (2.1) 
2 2 4. 

where M is the mass of the scalar field, and g is the Yukawa coupling. 

We split the field $ into its expectation value and the quantum fluctu- 
ations r\: 

$(x,t) =0(t)+77(x,t) , (2.2) 

with 

0(t) = ($(x,t)) = ^|l. (2.3) 



The scalar fluctuations have already been analyzed in |]33|, 53|. The 



equations for the system we consider here, with the back- reaction of the 
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fermion field, have been derived by |H] using the Schwinger-Keldysh formal- 
ism [^7], 28[ and the tadpole method fiS[| . We do not repeat it here. The 



equation of motion for the classical field are given by 

4>(t) + M 2 <P(t) + ^ 3 (t) + ^( V 2 ) + g^iP) = . (2.4) 

Here (ipip) and (if) are the expectation values of the quantum fluctuations 
of the fermions and the scalar field, respectively. They are related to CTP 
(closed-time-path) Green functions. They can be expressed by mode func- 
tions which satisfy the linearized equations of motion in the background field 
and initial conditions at some time to- in the following we choose to = 0. The 
scalar back-reaction via (rf) has been calculated previously by various groups 
within different approximation schemes, among them large N, Hartree, or the 
one-loop approximation. As we have explained above, here we are merely in- 
terested in the fermionic back-reaction and do not include the scalar one, 
except for some of the numerical examples in section 8. 
The fermion field ip satisfies the Dirac equation 

[idt-H(t)]i/>(t,x) = Q, (2.5) 

where the Hamiltonian Ti is given by 

H(t) = -iaV + m(t)(3 . (2.6) 

The term m(t) = g<fi(t) is the time dependent fermion mass. We expand the 
fermion field in terms of the spinor solutions of the Dirac equation 

X ) = E/ TTtLt [&P..EW*) + ^ PiS V- PlS (t)] e +ip ' x , (2.7) 

with the time independent creation and annhiliation operators whose mass 
is determined by the initial state. The creation and annihilation operators 
satisfy the usual anti-commutation relations 

{b P ,s,bl,A = 2E (27i) 3 6 3 (p-p')6 ss , , (2.8) 
{d P ,„ dj,,-'} = 2E (2rt) 3 5 3 (p-p , )5 ss/ . (2.9) 

For the positive and negative energy solutions we make the ansatz 

U p , s (t) = N [id t + H p (t)} f p (t) ( X ° ) (2.10) 
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and 





Xs 



V p , a (t) = N [id t + H- P (t)]g p (t) 

with the Fourier-transformed Hamiltonian 

H p (t) = ap + m{t)(3 . 
For the two-spinors Xs we use helicity eigenstates, i.e., 

V°X± = ±X± • 



(2.11) 



(2.12) 



(2.13) 



The mode functions f p and g p depend only on p = |p|; they obey the second 
order differential equations 



d 2 

— - irh(t) + p 2 + m 2 (t) 



dt 2 



+ im(t) +p 2 + m 2 (t) 



hit) 
9 P (t) 



0, 
. 



(2.14) 
(2.15) 



The initial state for the fermion field is usually specified as a vacuum or 
thermal equilibrium state obtained by fixing the classical field 0, and thereby 
the fermion mass, to some value 0o for t < 0. The spinor solutions are 
then identical with the usual free field solutions of the Dirac equation with 
constant mass m Q = m(0) = g<fi(0). Therefore, the mode functions, which 
would be plane waves for t < 0, satisfy the initial conditions 



/ p (0) = l , f p (0) = -iEo, 



(2.16) 
(2.17) 



g p (0) = l , g p (0) = iE o . 
For the spinors U and V we use the usual free field normalization conditions 



C/ P)S (0)[/ P , S (0) = -y p , s (0)Vp )S (0) = 2m 
C/ PiS (0)C/ p , s (0) = ^(0)^(0) = 2E Q ; 

we will also need the orthogonality relation 

^.(0)V-p^(0) = . 

i?o(p) denotes the mode energy in the initial state, 



El 



p 2 + m 2 



For the normalization constant we find 

iV = [£ + mon 1/2 . 



(2.18) 
(2.19) 



(2.20) 



(2.21) 



(2.22) 
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For t > Eqs. (g3g ), (gT| ), and (|2TT6| ) imply that 

f P (t) = g* p {t) . 



(2.23) 



Since the time evolution of the spinors U P)S {t) and V- PtS (t) is induced by 
the hermitean operator 7i p , their normalization and orthogonality relations 
(|2.18| ) and ( |2.20| ) are conserved. This implies a useful relation for the mode 
functions |yj 

\f p (t)\ 2 -im(t) [f p (t)f;(t) - f p (t)f;(t)] + [p 2 + m 2 (t)] \f p (t)\ 2 = 2E (E + m ) 

(2.24) 

which takes the role of the Wronskian. Using these mode functions (ip?p) can 
be calculated once the initial state is specified Q. If we use the Fock space 
vacuum defined by b PtS \0) = and d PtS \0) = we get 

d^p 



E 



(2n)Z2E V - p > 



d^p 



(2vr) 3 2£o 



2E, 



2p 2 



E + m 



\f P \ 



(2.25) 



If we use a thermal density matrix defined in terms of the Fock space states 
one obtains 



#P . ,/£o\L 2p 2 , f 
(27r)' :i 2_Eo V2T/ I So + mo 



(2.26) 



the integration measure in the momentum integrals is modified accordingly 
in the expressions for the energy-momentum tensor. 
We will denote (ipip) as the fluctuation integral 



m = mm)) 



(2.27) 



The fluctuation integral is divergent and has to be regularized and renormal- 
ized. This will be done in section 4. 



3 Energy momentum tensor and particle num- 
ber 

The energy momentum tensor of a Dirac field with mass m = m(t) is given 

^ 1 ~ 

= i> f-r^ d v +mg^ ip . (3.1) 

4 We use the Heisenberg picture, i.e., the field operators depend on time via the mode 
functions. 
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The expectation value of the energy momentum tensor taken with the initial 
density matrix is spatially homogeneous and has therefore the form = 
diag(£,"P). In addition to the contribution of the Dirac field it contains the 
classical contribution of the scalar field. The energy density £ of the quantum 
fluctuations is then obtained as 

£ a (t) = 08Kp) V)=E/ { 2n%E V ' pAt) m ' p) V ~ pAt) (3 - 2) 

= 2 / (2n%E ^ [E ° " m ° )] " ^) " 2Eom( ~ t )} ■ 

Using the equations of motion it is easy to see that the time derivative of the 
total energy density 

£ = £* + &= l -tf{t) + l -M 2 <p\t) + ^<p\t) + £a(t) (3.3) 
vanishes. The fluctuation pressure is given by 

Mt) = I(^ ) = ix:/^^F_ p „( t )7^-p„W (3.4) 

= 5/ (2T ^ 2Eo {[ g ° - m »] [< (a/; - i/;) - m«)i/„i 2 ]} , 

the total pressure is 

V(t) = <p 2 (t)-£ + V R . (3.5) 

Energy density and pressure are quartically divergent, their renormalization 
will be discussed in section 4 along with the renormalization of the fluctuation 
integral. 

In contrast to the fluctuation integral and the energy-momentum tensor, 
the definition of the particle number relies on the creation and annihilation 
operators. The number of particles with momentum p and helicity s is given 
generally via 

KAt) (blMKAt)) ■ (3.6) 

The definition of time-dependent creation and annihilation operators implies 
an interpretation f\. If we used the operators 6 PjS and the particle num- 
ber would remain equal to the initial one, zero for an initial vacuum state, 
or 1/ [exp(E /T) + 1] for a thermal one. While these operators refer to a 
decomposition of the field with respect to the exact mode functions f p (t), 
the concept of free particles implies plane wave mode functions. If we define 

5 For an extensive discussion, in the context of general relativity, see [Eol . 
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these modes to have mass m and Eq = \Jp 2 + Mq, i.e., if we use the modes 
corresponding to the initial state, we obtain 

b° p>s (t) = J d'xe-^^M^t) 

= Cl s (t)bp,s + D ptS (t)dl p>s (3.7) 

with 

C= 7 ^[^ + «p(*o)](^). (3.8) 
We need only the coefficient D p s which is given by 



p iE t 

D° p , s = 2^"^-p..(*) > (3.9) 
and the particle number becomes 

KAt) = \Dls(t)\ 2 - ( 3 - 10 ) 

In terms of the mode functions we obtain for the occupation number for one 
helicity eigenstate with momentum p 

K^) = [2E + i[f;(t)f p (t) -f;(t)f p (t)] -2(mx -mo)\f P W 

(3.11) 

This definition has been used in As it should be for fermions, the 

occupation number is strictly less or equal to 1; this is obvious from Eq. 
(|3.9|), since D ps is the scalar product of two complex vectors of unit norm. 
Integrating the occupation number over momentum one obtains the total 
particle number density 

W) = E/ j^KAt) ■ (3-i2) 

If we imagine the time evolution being stopped at the time t it seems more 
natural to use free quanta of mass mi = m(t) and 

b P ,M = J d 3 xe-^ + ^U p ]M^t) (3.13) 

with 

^«=7irW Bl+ ^ (tl)1 (o)- < 3 - 14 > 
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The coefficient D^ s now reads 

D *--=m u " v -* Atl) ' (3 - 15) 

and the occupation number for particles with mass m\ = m(t) becomes 

E - m 



4E E 1 (E 1 + m 1 

+i( El +m 1 )[f;(t)f p (t)-f;(t)f p (t)] 



2E (E + m ) 

(3.16) 



The total particle number J\f\(t) is again obtained by integrating J\fp S over 
momentum and summing over helicities. The particle number is divergent by 
power counting, the analysis of the divergent contributions of M\(t) shows, 
however, that it is finite and does not need counter terms (see also section 
4). For Ao(t) we find a linearly divergent contribution that vanishes in di- 
mensional regularization. 

As we have mentioned above, the definition of the particle number relies 
on an interpretation, which seems to be more straightforward if the particle 
number is computed for particles with the "final" mass mi = m(t). Of course, 
even if the classical scalar field relaxes to as t — > oo, the final state (taken in 
the Schrodinger picture) never becomes an ensemble of free particles. Such a 
state would be described by a density matrix which cannot arise in unitary 
evolution from a pure state |0). This is also true for the additional particles 
created in the case of a thermal initial state. 

Eq. fl2.24p may be used to recast the expressions for energy density, 
pressure and particle number into a different, sometimes advantageous form. 



4 Renormalization 

In order to develop the framework for renormalizing the one-loop equations, 
we write the equation of motion for the mode functions, Eq. (|2.14| ), in the 
form 

i] ' ' f P (t) = -V(t)f 9 (t), (4.1) 



with 



V(t) = m 2 (t) -ml -im{t) . (4.2) 

Using the initial conditions ( |2.16| ) this equation can be recast into the form 
of an integral equation: 



fv(t) 



-iE t 



E Jo 



dt'sm[E (t-t')}V(t')f p (t') . 



(4.3) 
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Using this integral equation, the mode functions may be expanded with re- 
spect to the potential V(t). We split off the zeroth order (plane wave) con- 
tribution and an oscillating phase factor by writing 

f p (t) = e-^ t [l + h p (t)}- (4-4) 

The functions h p (t) satisfy differential and integral equations derived from 
Eqs. ( |4.1| ) and (|4.3|), respectively. These functions are discussed in Appendix 
A. They may be decomposed as 



h p {t) = Y,h p n \ (4.5) 

n=l 

where h p n > is of n'th order in V(t); we define further the inclusive sums 

oo 

h p n) = E h p m) ■ (4.6) 



p 

m=n 



In terms of these functions and their expansion discussed in Appendix A the 
integrand of the fluctuation integral can be written as 

l-^)l/ P W! 2 = ^-(l-^){2Reh p (t) + \h p (t)\ 
-C/n / ran \ "^n / 



m V m 
m(t) m(t) m 3 (t) m(t)m 2 (0) 

~E7 ~ " 2{Eof + 2(E ) 3 

+^^cos(2E t)+K F (p,t). (4.7) 

The first terms on the right hand side lead to divergent or singular momentum 
integrals. The function Kp(t) can be considered being defined by this equa- 
tion. It behaves as (E Q )~ 4 and its momentum integral is finite. While Kp(t) 
is defined here as the difference between the original, numerically computed 
integrand and its leading contributions, an alternative expression, avoiding 
such a subtraction, is given in Appendix A. We decompose the fluctuation 
integral as 

Ht) = ^divW + ^sing(t) + ^fin(t) , (4-8) 

with 

r (ftp j m(t) rh(t) m 3 (t) m(t)m 2 (0)l 
div = l^\^~W^~W*f + 2{E f /' (4 ' 9) 



9 



The divergent part Tdn is a local polynomial in m(t) = g<p(t). It will be ab- 
sorbed by appropriate renormalization counter terms. The integral involving 
cos(2E t) is logarithmically singular at t — but finite otherwise. This con- 
tribution is obviously related to the initial conditions and will be discussed 
in section 5. 

The fluctuation parts of the energy momentum tensor can be analyzed in 
a similar way. The integrand of the energy density can be expanded as 



^) (/„/; - fpfi) - m(t) 



= -(Eo-mo) {l + 2Re/i p + \h p \ 2 + -^Im [h p (l + J - m(t) 



— —En 



m 2 (t) ml m 2 (t) m 4 (t) 
2E + 2E~ + 8{E ) 3 + 8{E y 



8(£ ) 3 4(E ) 



+ K E (p,t) . 



(4.12) 



Again K E (p, t) is defined by this equation and it behaves as (E ) 4 as E — > 
oo. There is no cosine term here and, therefore, no singular contribution. So 



£{t) = E& v {t) + £ &n (t) , 



(4.13) 



with 



div 



■'fin 



<fip 



-Eq- 



m 2 (t) ml m 2 (t) m 4 (t) 
~2Eo~ + 2E~ + 8(£ ) 3 + 8(£ ) 3 



m\ m 2 {t)m\ 
3 ~ 4(£ ) 3 

K E (p,t) . 



8(£ ) 3 
(ftp 



(4.14) 
(4.15) 



If the integrand of the fluctuation pressure is rewritten in terms of the func- 



tions hp it reads 



-(£ -m ) + (l + 2Re/i p + |/g 2 ) - ^Im [hp (l + h* p )] |-m(*) . 

(4.16) 

We expand again in V(t) in order to sort out the leading contributions. 
Finally, V(t) can be decomposed as 



V(t) = V div (t) + V sing (t) + Vsn(t) , 



(4.17) 



10 



with 

r d^p j 4 2m\ m o rn 2 {t) m{t)m{t) m 2 (t) 

div = Jj^f \ 3 + ^ + 6(W + 6GW" 12(£ ) 3 ~W 

- 2 w-n (418) 
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^ /• tft> m(t)m(0) , « „„x 

Pfl ° = I Wf Kp{r ' t] ■ (420) 

The integral over .fTp, which is defined by this decomposition, is finite. The 
divergent terms JF div , £ div , and Pdiv are proportional to local terms in <f>(t) and 
its derivatives. These can be absorbed in the usual way by introducing the ap- 
propriate counter terms into the Lagrangian and into the energy-momentum 
tensor. 

The counter terms in the Lagrangian are introduced as 

C c , = - \bM 2 <? - ^ . (4.21) 

The divergent parts of the fluctuation integral can be evaluated, e.g., using 
dimensional regularization. One finds 

m(t) m 3 (t) m(t)m 2 
•fdiv - -^L + —L + , (4.22) 

with the abbreviation 

L = - + ln^- 7 . (4.23) 
e ms 



As already found for the scalar fluctuations |33| , the dependence on the initial 
mass m can be absorbed into finite terms, AZ, AM 2 and AA. Applying an 
MS prescription, the infinite renormalizations become 

5Z = -£- 2 L, (4.24) 
8X = -e£- 2 L, (4.25) 

with 

2 , Ann 2 

L = 1 + ^ 1 ^r-l- (4-26) 
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There is no infinite mass renormalization counter term. Introducing the 
renormalization counter terms into the equation of motion, we obtain 



A + AA 
IT 



1 + AZ) 4> + [M 2 + AM 2 ) 4> + ^—<f + g{T^ + ^ sing ) = . (4.27) 



The coefficients of the finite terms left over after adding the renormalization 
counter terms to jF div are given by 

q 2 M 2 

AZ = ^ln— , (4.28) 
q 2 M 2 

AA = 6-^4 In— , 4.29 

All 1 TTIq 

AM 2 = (4.30) 

Since the bare fermion mass vanishes, we have introduced the scalar mass M 
as scale parameter. Obviously, the equation of motion is not yet acceptable 
in its present form, due to the singular term. 

The divergent parts of the energy give, after dimensional regularization 



™ 2 (t) T m A {t) m 4 (0) m 2 {t)m 2 , 
lo7r lb7r z 62-k z 8tt z 

The counter terms correspond to those in the Lagrangian, i.e., 



£ c , = + \bM 2 <? + (4.32) 

with the same coefficients as above. We need no infinite counter term for the 
zero point energy or cosmo logical constant. Adding the divergent part and 
the counter terms we are left with finite contributions 

£ div + E cX , = l -AZtf + l - AM 2 <p 2 + |^0 4 + AA , (4.33) 



with 



^ = -^4 ■ (4.34) 



32tt 2 

The divergent part of the pressure is given by 

_ m 2 {t) m(t)m(t) ml m 2 m 2 (t) 

Vdiv ~ l2^ L ° ~ 24^ Lo ~24^ + ~Yh^ ■ (435) 
In addition to the counter terms already introduced we have to add to the 



energy momentum tensor the "improvement" counter term |39 

5A(g^d a d a - d,d u )<p 2 (x) ; (4.36) 
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since (p depends only on t this term contributes only to the pressure, and we 
have 

V c . t = + 8A^f . (4.37) 

We choose 

M-^Ii (4.38) 

there is a finite remainder 



with 



q 2 M 2 

AA = -f-j In —r . (4.40) 



487r 2 m§ 



5 Removing the initial singularity 

We are now ready to discuss the terms proportional to cos(Eot) which turn 
out to be singular as t — > : 

= Ifja cos ( 2 ^) - ln(2-ot) as * - . (5.1) 

For the case of scalar fields in Minkowski space the fluctuation integral has 
only a logarithmic cusp oc tint at t = 0, as observed by Ringwald Q; the 
energy is finite, and the pressure behaves as Int. In FRW cosmology the 
energy of scalar fluctuations is logarithmically singular whereas the pressure 
behaves as 1/t. So the Friedmann equations become singular. Problems 
with the initial conditions in FRW cosmology have also been noted in 0, [L0[| 
when using comoving time and the associated vacuum state; the problems 
disappear if conformal time is used. The two vacuum states are related by a 
Bogoliubov transformation. 

For fermionic fluctuations we find that already the fluctuation integral is 
divergent, so that the numerical code cannot be started even in Minkowski 
space. We have shown recently, for the case of scalar fields, that such "Stueck- 
elberg singularities" |35| can be removed by a Bogoliubov transformation of 



the initial state, which was constructed explicitly. 

Within the Fock space based on the "initial vacuum" state |0) which is 
annihilated by the operators b PyS and d PjS we define a more general initial 
state by requiring that 

"&p,.-Pp,.d!pJ|6) = 0. (5.2) 
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The Bogoliubov transformation from |0) to this state is given in Appendix 
B. If the fluctuation integral, the energy, and the pressure are computed by 
taking the trace with respect to this state, we just have to replace in the 
defining equations (|2.25| ), (|3~2j ) and (|3.4j ) the functions U PjS (t) by 



V_ PiS (t) => cos(/3 PiS )y_ PiS (t) + sin(P P;S )U PiS (t) 



(5.3) 



For the particle number, the substitution is done in the Bogoliubov coeffi- 
cients, Eqs. ( p.9|) and ( |3.15| ). The angle (3 P)S is related to p P;S via 



p PiS = tan(/3 PjS 



(5.4) 



If the expectation value of ipip is taken in the Bogoliubov-rotated initial state 
the fluctuation integral becomes 

+ \U P:S (t)V^ s (t) + V P:S (t)U. PiS (t)]smp P:S cosp PiS } . (5.5) 
Rewriting this expression in terms of the mode functions we find 



E 



cPp 



+ sin2/3 Pi 



(2tt) 3 2£ 
sp 



cos 2/5 Pv , 



2E ° ~~ ~E^~7 I^pI 2 

E + m 



E + m 



Imd t f 2 -2m(t)Ref 2 



(5.6) 



Using the perturbative expansion of the mode functions given in Appendix 
A this integral takes the form 



E 

s 

m(0) 



cos 2/3, 



p.* 



m(t) m(t) m 3 (t) m(t)m 2 (0) 



E 4(£ ) 3 2(E ) 



2(E^ 



+ sin2/3 p 



3 cos{2E t) + K F {p,t) 
sp 



-2E cos 2E t + L F (p,t)} 



•2E (E + m ) 

The functions K F (p,t) has been defined above, L F (p,t) is defined by 



(5.7) 



L F (p,t) = 21m 



2m(t)Re/ 2 . (5.8) 



~e- 2iEot ((1 + hp% - iE (2h p + fi: 

Obviously, one gets rid of the term proportional to cos(2E t) by requiring 

m(0)(E + m ) 



tan 2j3 p <s 



8sp(E ) 



(5.9) 
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Thereby, the Bogoliubov transformation is explicitly specified. We notice 
that the helicity dependence of the Bogoliubov transformation cancels in Eq. 
(|5.7|). Using the asymptotic behavior 



p— >oo 



m(0) 
8sp{E f 



and therefore 



cos(2/3 p , s ) - 1 =2sin 2 (/3 PiS ) P ~°° 



m 2 (0) 



sin(2/3 Pi . 



p^oo 



m(0) 



(5.10) 



(5.11) 



(5.12) 



4sp(£ ) 2 ' 

it is easy to convince oneself that this Bogoliubov transformation does not 
interfere with the analysis of the divergent parts and, therefore, with the 
renormalization discussed in the previous section. So T is rendered finite by 
adding the counter terms defined in the previous section. In the renormalized 
equation of motion ( |4.27| ) we just have to replace ^ S i ng (t) + TfaS$) with ^g n (t) 
which is given explicitly by 



^fin(t) 



-2 sin 2 (3, 



p..s 



m(t) m(t) m 3 (t) m(t)m 2 (0) 



+ cos 2j3 PtS K F (p, t) + sin 2/3 Pi , 



4(^)3 
sp 



2E (E + m 



2(£ ) 3 
-L F (p,t) 



+ 



2{E ) 



(5.13) 



The renormalization of the energy density proceeds in an analogous way. The 
fluctuation energy in the Bogoliubov transformed state is 



E 



d^p 



{cos 2 /5 p ., V^ s {t) (PH P )V- P , s (t) 



(2tt)32£ 

+ sin 2 (3 p , s U p , s (t)(f3H p )U^ s (t) (5.14) 
+ sin/5 p , s cos(3 p , s [u PyS (t) (/3H P ) V„ PjS (t) + V Pt ,(t) (/3H P ) U- Pt8 (t)] } . 

We again insert the expansion of the mode functions to obtain 
d^p 



E 



(2vr) 3 2£ 



{cos2/5 p , s [i [E - m(0)] (fj; - f p f;) - 2E m(t) 



+ sin2/3 Pv 



sp 

E + m ' 



-Re 



ft + E 2 f? 



(5.15) 
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After adding the counter terms defined in the previous section the finite part 
of the fluctuaion energy becomes 



-2sin 2 /? p „ 



m 2 (t) ^rn 2 m 2 (t) 



2Eq 2Eq 8(-£/q 



+ 



m 2 (t)ml 
+ cos 2/3 P!S K E (p, t) + sin 2/3 Pj , 



(5.16) 



2£ (£ + m ) 



with 

L E (p,t) = Re {e- 2 ^ 



• 2 



h p - 2iE h p (1 + h p ) + (m 2 (t) - m 2 )(l + h P Y 



}■ 

(5.17) 



Finally, we consider the fluctuation pressure, taken in the new vacuum state. 
It reads 



(fip 
(2^ 



cos 2/? P)S 



^ (fj; - f P f;) (Eo - m ) 



--m(t) {E - m ) \f p \ 2 - 2m{t)E 



(5.18) 



+ sin2/? p 



sp 



-Re 



\fl-^mmfl+{\v 2 + \rn\t))f 2 p ]) . 



p ' s E + m " 

The finite parts becomes, after adding the counter terms, 



-2sin 2 /3 Pv 



(2tt) 3 k 

m(t)m(t) m 2 (t) m 2 {t)rn^ 
~ 12(£ ) 3 ZE 6(£ ) 3 

+ cos 2f3 P:S K P (p, t) + sin 2(3 P 



4 2m 2 , m 
~3 0+ 3E~ + 6(Eo) 3 ' 6(E ) 3 



(5.19) 



2£ (E + ™o) 



L P (p,t) 



with 



L P (p, t) = ^Re {4/ 2 - 2ime- 2iEot [h p {l + h p ) - iE (2h p + h 2 p ) 



+ 



A P 2 + 2m 2 (t)]f 2 } 



(5.20) 



The Bogoliubov transformation has removed the singular term in the pressure 
as well. 
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6 Extension to FRW spacetime 



Now that we have set all basic equations and performed renormalization the 
extension to FRW spacetime is straightforward. 

We consider the Friedmann-Robertson Walker (FRW) metric with cur- 
vature parameter k — 0, i.e., a spatially isotropic and flat spacetime. We 
will treat the the quantum fields and the cosmological background self- 
consistently. That is, the scale parameter a(t) is obtained dynamically from 
the quantum fields. 

The line-element of a flat FRW universe is given by 

ds 2 = dt 2 -a 2 (t)df 2 . (6.1) 
The time evolution of the a(t) is governed by Einstein's field equation 

Gfiu + aH$ + /3Hj£ + Ag, u = -k(T^) (6.2) 
with k = 8ttG. The Einstein curvature tensor G^ v is given by 

Gfw — R^u — ^9 V-^ • (6-3) 

The Ricci tensor and the Ricci scalar are defined as 

Rfw = R^, u x > (6-4) 

R = g^Rp,, (6.5) 

where 

R a/3~, = ^7 r a/3 _ + T^T^ - T^T^ . (6.6) 

The terms H$ and H$ arise if terms proportional to R 2 and R^R^u are 
included into the Hilbert-Einstein action. If space-time is conformally flat, 
these terms are related by 

H$> = I HjS , (6.7) 



so that we can set (3 = in (|6.2|) without loss of generality [f|(J. We also 
replace H^J by H^ u in the following. 

These terms are usually not considered in standard cosmology. They 
are included here, as well as the cosmological constant term, only for the 
purpose of renormalization; they will absorb divergences of the energy mo- 
mentum tensor. So in principle they should appear on the right hand side 
as counterterms; they are related to the coefficients of these counterterms by 
A = k5A and a = n5a. 
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As usual we can reduce the Einstein field equations to an equation for the 
time-time component and one for the trace of G^ u , the Friedmann equations 



Gtt + aH tt + A 
G£ + aH» + 4A 



(6.8) 
(6.9) 



-3H 2 (t) , 
-R(t) , 

-6 (H{t)R(t)+H\t)R{t)-^R\t) 



For the line element (|6.1|) the various terms take the form [|40 

G tt {t) 

H tt (t) 

H£(t) = -6(R(t) + 3H(t)R(t) 
with the curvature scalar 

R(t) = 6 [H(t) + 2H 2 (t) 
and the Hubble expansion rate 

a(t) 



H(t) 



a(t) 



(6.10) 



(6.11) 



(6.12) 



The Dirac equation in FRW spacetime (see, e.g., [[40], [11]]) is given by 



id t +4^rl + 4tt«V - g<f>(t)y \ ^(t, x) = 
2a[t) a[t) J 



(6.13) 



It proves convenient to introduce conformal time and scales. The conformal 
factors for the scalar field and the fermion field are: 



?/>(£, x) = a~ 3/2 (t)^(r,x) 
0(t,x) = a(t) _1 0(r,x) , 



(6.14) 
(6.15) 



with x = a(t)x and dt = a(t)dr. In conformal time, and using these redefi- 
nitions of the fields, the Dirac equation simplifies to 



{id T + iaV - g4>{r)(3} ^(r, x) = 
If we introduce, as in section 2, the Dirac Hamiltonian 

H(t) = -iaV + g{p{r)P . 



(6.16) 



(6.17) 
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this equation takes the standard form 



id T -H(r) ^(r,x) = . 



(6.18) 



So the formalism of quantization can be taken over from section 2, if the other 
functions and operators are understood as being rescaled quantities as well. 
The rescalings for the other quantities of interest are, omitting arguments, 
super- and subscripts: 



x = ax 
b = ab 
E = a^Eo 
U = a-^U 

n = a- x H 



v 
d 



a x p 
ad , 



m = a 1 fh 



(6.19) 



V = 
N 



arWV , 
= a V 2 N 



The potential V(t) , Eq. ( [4.2| ), becomes the analogous expression V(r) with 
m(t) replaced by m(r). This means that the entire perturbative expansion 
and the analysis of divergences proceeds in perfect analogy to the Minkowski 
space analysis. The metric does not appear in this formalism, and therefore 
there are no divergences related to the metric. So the Dirac field does not 
contribute to the wave function renormalization of the gravitational field SZ g 
(or, equivalently, the renormalization of Newtons constant) and to the terms 

In the equation of motion the fluctuation integral scales as a -3 and so do 
the kinetic and the A</> 3 term. Therefore, the divergent parts of the fluctuation 
integral are absorbed by exactly the same counter terms 5Z(f> and 5\<j) 3 /6 as in 
Minkowski space. The same holds true for the finite remainders proportional 
to AZ and AA. However, the finite mass renormalization AM 2 is now 
replaced by a~ 3 AM 2 (p while the genuine mass term scales as a _1 M 2 0. 

The renormalized equation of motion for the scalar field <fi takes therefore 
the form 



d 2 ~ 

[l + AZ)— (f) + a 2 
dr z 



M'+(Z--)R 



+ AM 2 <b + 



A • A A - ., ~ 

^ 3 + ^fin = 0. 



6 



(6.20) 

The energy-momentum tensor generated by the fermionic fluctuations 
scales exactly as a~ 4 . 



7 The linearized equations of motion 

A simple intuitive approach to the interplay between the classical Higgs 
field oscillating with a frequency of the order of the Higgs mass M and 
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the fermionic fluctuations is to treat the system as a Higgs field decaying at 
rest into the fermions. For large amplitudes of the Higgs field this picture 
is certainly inadequate. For small amplitudes, e.g., at the end of inflation, 
the linearization of the equations of motion leads indeed to an approxima- 
tion which supports the simple decay picture. We will compare the exact 
equations and the lowest approximation in our numerical examples. Here we 
analyze the behavior of the system analytically along the lines in pi]. 

If we retain only the 0(g 2 ) part of the fluctuation integral, we find a di- 
vergent part which is cancelled by the wave function renormalization counter 
term, a finite part remains and the linearized equation of motion of the clas- 
sical field takes the form 



(1 + AZUCt) + M 2 0(t) + / dt'ZJt - t') (t') + £ r (tW0) = . (7.21) 

Jo 

The self-energy insertion is subtracted and is given explicitly by 

£ r (M') = -2g 2 J ^L-L cos (2E (t-t')) . (7.22) 
We define the Laplace transform of the condensate <fi(t) via 



dte~ st (j)(t) 



(7.23) 



the inverse transformation being given by 



1 rioo+c 
Z7TZ J — 2oo+c 



(7.24) 



The "Bromwich" -contour of the latter integral runs parallel to the imaginary 
axis. The constant c has to be chosen in such a way that ip(s) is analytic 
for Res > c. In our application ip(s) will have cuts along the imaginary axis 
and poles in the half-plane Res < 0. 

For the Laplace transform ip(s) the equation of motion reads 



AZ) s 2 ip(s) - s<f)(0) - 0(0) 



(7.25) 



+AfV(s) 



-0(0) - s0(O) - s 2 0(O) + s 3 tp(s) + S r (s)0(O) = 



where S(s) is the Laplace transform of the self-energy kernel. The equation 
can be solved readily with the result 



#0 = 



0(0) + s0(O)l [l + AZ + st r (s) 



(1 + AZ)s 2 + M 2 + s 3 t r (s) 



(7.26) 
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The singularities of the right hand side in the complex s plane are given, 
on the one hand, by the singularities of S r (s), and by possible zeros of the 
denominator, whose locations have to be determined. Since 4>(t) cannot 
contain contributions that increase exponentially, these poles will have to lie 
on the left part of the complex plane. The dominant behavior at large t will 
be governed, therefore, by the singularities of E r (s) on the imaginary axis. 

The self energy kernel has been defined in Eq. (|7.22|) . Its Laplace trans- 
form is given by 

^w-^l^rar (7 ' 27) 

It has cuts along the imaginary s-axis running from s = 4imo to s = zoo, 
and from s = —4irriQ to s = —ioo. We introduce the frequency u = —is. 
Expressed in the variable u the cuts run from 4mo to oo and from — oo to 
— 4m . The discontinuity across the cut along the positive imaginary axis is 
defined as 

p(s) = disc£ r (s) = Ti r (iu + e) — S r .(zc<j — e) . (7.28) 

One finds 



, d 3 p 1 2ncu e . , g 2 JuJ 2 -Aml 

P 00 = -V - 2£o) = -7Z ■ (7-29) 



(2vr) 3 AEl ou + 2E K u; ' 4vr 



A more relevant quantity for ip(s) is the discontinuity of the denominator of 
Eq. O): 

,2 



«7(sJ 



discs ^(s) 



i^^/to 2 - Ami ■ (7-30) 



The discontinuity of the denominator is purely imaginary, as to be expected. 
If E is small, we can expand the denominator around its zero at s ~ iM. The 
contribution of this pole is related to the decay of the condensate particles 
of mass M into fermions of mass m$. We neglect the real part of £ as it just 
shifts the value of M. Indeed it should be zero at s = ±z'M if the condensate 
field is renormalized on shell. We write 

s 2 + (M-iT/2) 2 ~ -(uj 2 - M 2 + iMY) ~ s 2 + M 2 - f 7 (M)/2 . (7.31) 

It follows: 
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, M l{M) = ^M 2 -Am 2 . (7.32) 

This is almost, but not quite, what one would expect for fermions. Evaluating 
the width of <f> — > ff in the standard way one finds 

r^ /7 =^(M 2 -4^) 3/2 . (7.33) 
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A factor p is typical for an s-wave, the factor p 3 in the correct decay formula 
arises from the Dirac traces. While at first glance the difference is surprising, 
one may recognize that the factor, by which the two expressions differ, is 
given by 1 — Arn^/M 2 . The fermion mass is, however, of order g 2 , so that the 
difference is higher order and our approximation to the full one-loop result 
is just lowest order in g 2 only. 

The late-time behavior of <f>(t) is determined by the strongest singularities 
of its Laplace transform. These are, on the one hand, the poles at s — 
±zM — r/2 and, on the other hand, the branch cuts on the imaginary axis. 
The contribution of the poles, which are actually in the second Riemann 



sheet, has been analyzed carefully in plfl . If the poles have residue R they 
contribute 

~ —^iR (e m - v ' 2t + e~ lM - r / 2t ) = Re- Vt ' 2 2 cos(Mt) . (7.34) 

Approximately, R ~ 0(O)/2. While this contribution decreases exponentially, 
the singularities on the imaginary axis yield a power behavior 

<f)(t) = Rt~ a cos(ftt + if) , (7.35) 

where the power-like decrease as t~ a is related to the order of the branch 
point — 1 + a. Since the branch point is of the square root type we have 
a = 3/2. Our treatment differs slightly, in terms of order g 2 , from that of 
where a = 5/2. 

In our numerical computation we just find the exponential damping (|7.34j) 
with a value of T which agrees with the theoretical expectation ( |7.32| ). The 



power behavior (|7.35|) is apparently suppressed due to a small coefficient. 



8 Some numerical results 

We have implemented numerically the formalism developed in the previous 
sections. We will discuss in this section some results of our numerical simu- 
lations. We will pay special attention to the phenomenon of Pauli blocking 



invoked in |31]]. We have already seen that the occupation number can- 
not exceed unity on account of the unitary evolution of the mode functions 
U Pi s(t),Vp >s (t). So an unlimited parametric resonance cannot develop. A 
priori this should not limit the production of particles as the available phase 
space is large. However, by a phenomenon similar to parametric resonance, 
the production of particles turns out to be concentrated within a very small 
band and it is only in this resonant region where Pauli blocking can be ef- 
fective. One should keep in mind, however, that even for the bosonic case 
particle production shuts off in the one-loop approximation pU|. 
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We expect Pauli blocking to be especially effective if the initial amplitude 
of the inflaton field is large. A typical case is displayed in Figs. 1-5. It corre- 
sponds to the parameters M — 1, A = 1, g — 2, and 0(0) = 0.6. 0(0) is taken 
to be zero in all examples. We show the behavior of the inflaton amplitude 
in Fig. 1, the conserved total energy and its classical and fluctuation parts in 
Fig. 2 and the pressure in Fig. 3. The time dependence of the particle num- 
ber is displayed in Fig. 4, using both definitions: Ao(t), referring to quanta 
of mass rriQ and J\fi(t) referring to quanta of mass m(t). The latter one is 
seen to behave more smoothly. The momentum spectrum of the occupation 
number J\fp S (t) varies strongly with time. We display therefore, in Fig. 5, 
the envelope obtained by selecting the maximal occcupation number reached 
at fixed momentum, as a function of Eq. The structure of this envelope 
shows resonance-like enhancements at threshold Eq = 1.2, and at Eq = 1.9. 
The maximal occupation number is reached only at these two values. On 
the one hand there is no Pauli blocking in the sense that all levels would be 
maximally occupied, on the other hand the unitary evolution does not allow 
for a parametric resonance with high occupation numbers. Clearly, for this 
parameter set the classical field is not able to efficiently transfer energy into 
the fermion fluctuations and its amplitude stays essentially constant. These 
results are similar to those obtained in |31] . 



The situation is simple in the case of very small excitations and moderate 
couplings. In this case the scalar field can decay into fermion-antifermion 
pairs. An exponential decrease is found for the exact quantum evolution and 
for the linearized equations of motion. This is displayed in Figs. 6-10. Fig. 6 
shows the exact evolution of <ft(t) for the parameter set M — 1, g — 2, A = 1, 
and 0(0) = 0.01. The amplitude is seen to decrease exponentially, the decay 
rate being given approximately by T ~ g 2 /87T. In the same Figure we also plot 
the solution of the linearized equation of motion for which T is exactly equal 
to g 2 /8ir. The energy is transferred completely to the quantum fluctuations. 
The pressure, plotted in Fig. 7, becomes asymptotically equal to one third 
of the total energy density £ = 7.5 • 10" 5 . The quantum ensemble created 
is ultrarelativistic, as to be expected for massless quanta. The momentum 
distribution of the occupation number is shown in Fig. 8, it is characterized 
again by a resonance-like band. The total particle number is plotted in Fig. 
9. 

So far the results correspond to the expectations. The situation is however 
not as transparent. For intermediate initial amplitudes 0(0) the relaxation 
can shut off even for case 2m(0) < M where the scalar field can decay. An 
example is given in Fig. 10, with the parameters M = 1, g = 1, A = 0, 
and 0(0) = 0.1. The full evolution stagnates, the linearized equations of 
motion show the expected exponential decrease. On the other hand, even for 
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large initial amplitudes the transfer of energy can be as efficient as for scalar 
fields. An example is displayed in Figs. 11-13. The amplitude 4>(t), plotted 
in Fig. 11, decreases by roughly a factor 3. The energy transfer is evident 
from Fig. 12. The spectrum of the occupation number, presented in Fig. 13, 
shows that the resonance type band is occupied more strongly than in the 
first example, which had a smaller initial amplitude, and otherwise the same 
parameters. Obviously, the concept of Pauli-blo eking is too simple to describe 
the situation in an adequate way. If one aims at a better understanding 
of the quantum evolution, analytical methods should be developed. The 
simplest approach could be an analysis of the differential equation of the 
mode functions for a given oscillating classical field, analogous to the analysis 
of the Mathieu or Lame equations. In order to illustrate the behavior for a 
given oscillating field we have solved numerically Eq. ( 2.14j) with m(t) = 
mo cos t for various values of mo and momentum p. We plot in Fig. 14 
the envelope of the occupation number Af env , i.e., the maximal occupation 
reached at fixed function of E (p). The structure resembles the 

envelopes plotted in Figs. 5 and 8. 

Finally, we shall present an example where we include quantum fluctua- 
tions of both the fermion field and the scalar field itself. In Fig. 15 we display 
the behavior of the amplitude for this combined system, the parameters are 
M — l,g = 2, A = 4, and 0(0) = 1. The field is seen to relax efficiently. 
On the contrary, if the scalar fluctuations are not included, the relaxation in- 
duced by the fermionic fluctuations is small, as seen in Fig. 16. Fig. 17 shows 
the relaxation for the case that only the scalar fluctuations are included, it 
is seen to start later and to be less efficient than for the combined system. 
The growth of the fermionic and bosonic energy density is plotted in Fig. 18. 
The fermionic energy density is smaller but rises earlier. Its asymptotic value 
is only by roughly 20% higher than in the purely fermionic evolution. The 
fermion fluctuations seemingly act here as a kind of catalysator supporting 
the development of the bosonic quantum fluctuations. 



9 Conclusion 

We have developed the quantum field theory of the out-of-equilibrium evolu- 
tion of fermionic quantum fluctuations driven by a scalar field. The quantum 
back-reaction has been taken into account in one-loop approximation. We 
have formulated the renormalization of the equations of motion and of the 
energy-momentum tensor in a covariant form and independent of the ini- 
tial conditions. A restriction of suitable initial conditions for the fermionic 
quantum system, as required by the removal of initial singularities, has been 
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obtained by selecting a Fock space built on a Bogoliubov-transformed vac- 
uum state. Furthermore, we have formulated the renormalized equations for 
the case of a spatially flat FRW metric. 

We have numerically implemented the evolution equations and we have 
presented some examples for the evolution of the quantum system. If the 
initial amplitude of the scalar field is very small the system evolves as pre- 
dicted for the linearized equations of motion, formulated in section 7. It can 
be described as a decay of the scalar field into fermion-antifermion pairs. If 
the initial amplitude is larger, the evolution depends on the way in which a 
kind of resonance band at low momenta is situated kinematically and how 
it is filled. In some cases the fermions are indeed ineffective in damping 
the oscillation of the classical field, in others the relaxation develops as in 
the bosonic case and a considerable part of the energy is transferred to the 
quantum fluctuations. Analytical studies should help to clarify the features 
observed for large-amplitude oscillations. An example where bosonic as well 
as fermionic fluctuations are included shows an interesting interplay where 
the fermions catalyze the development of bosonic fluctuations. 

We think that these results show that nonequilibrium systems with fermionic 
fluctuations show more interesting features and may play a more interesting 
role in cosmology than previously assumed. 
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A Perturbative expansion of the mode func- 
tions 

We have introduced in sections 2 and 3 the mode functions f p (t) and h p (t) 
which are related via Eq. (|4.4j ). In this Appendix we will analyse the per- 
turbation expansion and ultraviolet behavior of the functions h p (t). These 
mode functions satisfy the differential equation 

hp - 2iE h p = -V(t) [1 + hp] , (A.l) 

with the initial conditions h p (0) = h p (0) = 0. We expand h p with respect to 
orders in V(t) by writing 

h p = hg ) +hP+h® + ...., (A.2) 

where h p n \t) is of n'th order in V(t). h p n ' denotes the sum over all orders 
beginning with the n'th one. 



^ B) = E^ B) . ( A - 3 ) 

l=n 

so that 

h p = = + hf . (A.4) 

The integral equation for the function h p can be derived in a straightforward 
way from the differential equation satisfied by the functions f p ; it reads 

kp = 2Eo lo dt ' i*"*™ " ^ V{t ' ] [1 + kp{t ' )] ■ (A - 5) 

Using this integral equation we can obtain the functions h p n \t) by iteration 
33| . h^ is given by 



/#) = [ dt' ( e 2iE ^ - l) V{t') . (A.6) 



2E Jo 

Using integrations by parts this function can be analyzed with respect to 
orders in E via 



+ (*?'). < A - 7 > 
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with 



2Ec 



n+1 ,. t 



dt , e 2iE {t-t')y{n)^ 



(A. 



Here V">(t) denotes the Zth derivative of V(t); the subscript n indicates that 
the expression in parentheses has been reduced to negative powers of E equal 
or higher than n. For energy density and pressure we need the expansion of 
h& (t) as well. From Eq. ( |A.8|) and the relation 



= 2iE h ( p - f dt'V{t') 



(A.9) 



we find 



E 

1=0 



\2Er 



l+l 



V (l \t) 



o 2iE t 



V {1 \0) 



—i \ n+1 rt 



(■21 



dt f e 2iE (t-t')y(n+l)^ 



(A.10) 



Similar expressions hold for the higher and hy\ 

In the numerical implementation the functions \ can be obtained 

as the Fourier transform of the ra'th derivative of V(t). Its computation need 
just one update per time step. Alternatively one may construct differential 
equations satisfied by these functions, e.g., 



2iE (/£>) 



2E n 



-V . 



(A.n; 



In order to isolate the divergent terms in the fluctuation integrals of the 
equation of motion and of the energy momentum tensor we need an expansion 
of the mode function up to order O \{Eq)~ 3 ] and O [(E'o) -4 ], respectively. For 
this reason we will give relevant expansions of the h p in the following. 

The expansion of h^' up to O [(E )~ 3 } gives 



2E Jo V ' 4(£, 



V(t) iV{t) iV(Q) 



+ 



8(£ ) 3 8(E 



S(E ) 3 Jo 



dt'e 2iE °^V(t) 



(A.12) 



For h!p we obtain 



1 



8(£ ) 2 

i 



dt' V(t') 



i 2 



+ 



8(S ) 



tdt'V 2 ^') 
Jo 



8(E ) 3 



V(t) j dtV(f) + (hW). 
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(A.13) 



where - includes all terms of h^> that have at least four negative powers 
of E . It satisfies the differential equation 



V(t) fdt'V{t') . 



(A.14) 



Finally, via integration by parts, takes the form 



^(3) = I 

P 48(£ ) 3 



£dt'V(t')\ 3 +(h^)_, (A.15) 



where the last term satifies the differential equation 



p Ji 



-2iE (hf). 



-v(t) 

lV 2 (t) rt 



iV(t) 



16(£ ) 3 Jo 
dt'V(t') . 



dt'Vit 1 ) 



(A.16) 



8(£ ) 3 

For the integrands of the fluctuation integral we need repeatedly the two 
expressions: 

JxGM) = 2Rehf + \hf\ 2 (A.17) 

and 

J 2 (p,t) = 2Im/i« + 2\mh^*h^ . (A.18) 
Using the expansion up to O [(E )~ 3 } we obtain 



Xi(p,t) = — dt'lmV(t')-^U + 



2(E y ^ 2(E y 



fdt'lmVit') 
Jo 



E Jo 

lmV(t) lmV(0) . „ . ReV(0) . „ , 
; - cos (2£ t) + \ J sin (2ff t) 



4(E )3 4(E )3 

1 rt 



A(E y 



4(£ ) 3 



/ dt'lmV 2 (t') + — — [ / dt'lmVU')] (A. 19) 



-^KeV(t)J*dt>IrnV(t>) + 2Re (hf)_+ , 
where we have introduced the terms of O [(-E )" 4 ] as 



2Re (hf)_=2Re (/£> + hf> + 



_ + 2Reh^ 



(A.20) 
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and 



J2> 
"P ' J a 



i 1 )|%+2Re(^)^)_+|^)| 2 + 2Re(^)^) . (A.21) 



Inserting the real and imaginary parts of the potential V(t) we finally have 
the result 



Zi(p,*) = 



-k [m{t) - mm -jh? 



m(t)m(0) -m 2 (0) 



+ 



1 



En 



\m z (t) + m 3 (0) - ^m(t)m 2 (0) 

-jm(0) cos (2-E7 f) + -m(0)m sin (2£ t) 
4 2 



cos(2£ t) 



+ 2Re (h^ + (\h^\% , 
Re V(t) lmV(t) hnV(0) 



(A.22) 



+ 



2(E y (2E Q ) 



ReV(t) , „ , ReK(t) ReV(O) , ^ , 

+ w sm{2Eot) + w - w cos{2Eot) 
+^ w sin(2Eot) + w Re m fo dt ' Im 

3[ReF(Q] 2 _ [lmV(t)] 2 _ lmV(t) r* 
+ 4(£ ) 3 4(£ ) 3 2(E )3 io im V[t } 

ReV(t) 



lmV(0) . 

sin 



Cdt'lmVit') 
Jo 



2(£ ) 3 

(2E Q t) fdt'ReVit') 

Jo 



+ Re Z^°) cos(2ff t) / dt'ReV(t') 



o 



2(E f_ 
+2Im (/ i W). + 2Im (h^h^ . 

For completeness we also give the expansion up to O [(E ) ~ 4 ]: 



(A.23) 



Ji(p,t) = ±r fdt'lmVit') 
£j{\ Jo 



ReV(t) 1 



2(£ ) 2 2(E )< 



dt' lmV(t') 



ImVCt) lmV(0) /orn , ReV(0) , , 
; cos (2£ t) + v / sin (2£ i) 



4(E )3 4(E )3 
1 



4(£ ) J 



I dt/M2(t ' ) + w[l dt ' M(t,) 



4(£ ) 3 ^o v ' 6(£ ) 3 

*B*v(t) /Vim no + 3[Re ;g )]2 - G^l! 

2(£ ) 3 W Vo 1 ; 8(£ ) 4 4(£ ) 4 
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]mV(t') 



4(E )4 
1 



f dt'lmV(t') - —l—r fdt'lmVit') f'dt'lmV 2 ^') 
Jo 4(_£/ ) 4 Jo Jo 



4(E ) 



ReV(t) 



dt' lmV(t r ) 



i 2 



Re 1/(0) 



4(£ ) 



cos{2E t) f dt' Re Wt') + 

JO 



24(£ ) 4 
Re frft) 

w 



4(£ ) 

imy(o) 



-ImV(0)sm(2E t) / (it Re V{t') - - \ J cos(2E t) 



o 



8(£ ) 



^sin(2E t) + 2Re (/,«)_+ (\h^\% , (A.24) 



where 



2Re (/*«)_ = 2Re (/£> + + /if + /if )_ + 2Re hf . (A.25) 

In terms of the notation introduced above the alternative representation for 
the functions Kp, Kp, and Kp become now 



K F (p,t) 
K E {p,t) 

K P (p,t) 



-E [l^p, + [J 2 (p, t)] 5 + m [Ji(p, f)] 5 
_mo 

4 
3 

m(t)mo 



(A.26) 
(A.27) 



[Za(p,*)]g , 

-_£J D [J 1 (p,t)] F + ^mo-^m(t)) [J, ( /,/)]- 



(A.28) 



4 4?7ln 
^ ^(p, f)fe + - [J 2 (p, - ^ [J 2 (p, • 



Unlike the case of scalar fluctuations [33], the evaluation of these expressions 
has required extensive algebra, and the numerical implementation requires 
further efforts. This reflects the fact, that fermion loops are divergent up 
graphs with four external lines. The definition of Kp, Kp, and Kp given in 
section 4 is much easier to handle, it involves, however, potentially dangerous 
differences between expressions which are computed numerically and leading 
perturbative terms. We found in our numerical computations that the simple 
subtraction was tolerable; this is due, here, to the fact that most integrals 
are dominated by the low momentum region. The analysis presented in this 
Appendix was necessary in any case, however, in order to find the divergent 
contributions. 



30 



B Bogoliubov transformation 

Here we recall the basic formulae for the Bogoliubov transformation of spin 
1/2 fields (see e.g. p2|). A general canonical transformation compatible 
with the anticommutation relations is given, up to a further trivial phase 
transformation, by 

V s = cos(/3p, s )6 PiS + sin(/3 PiS ) exp i,5p ' s rf f _ PiS (B.l) 
dt p>s = -sin(p P!S )e- i5 »'4 p>s + cos((3 p>s )d 1 _ P)S . (B.2) 

The state annihilated by the new operators b p>s and d P:S is given in discrete 
notation by 

|6> = II{cos(/3 p , s ) + — |— sin(/? p>s )e-^^ t _ p , s 6 PiS }|0) . (B.3) 

For the expectation values of the various bilinear products of interest here 
one finds, in continuum notation, 

(6|& P|S W|6} = (27r) 3 2E ^5 3 (p-p / )sin 2 /5 PiS , (B.4) 

(6|d_ P X P ', s '|6> = (27T) 3 2E ^5 3 (p-p / )cos 2 /5 PiS , (B.5) 

(0|gL p , s & pV |0) = (27r) 3 2E ^5 3 (p-p , )e i5p - s cos/3 PiS sin/3 PiS , (B.6) 

(6|fc P|S ^y|0) = (2vr) 3 2 J E;o^ ( 5 3 (p-p / )e"^cos/3 p , s sin/? PiS ,(B.7) 

so that the fluctuation integral becomes 

(2ir) 3 2E { Up > sUp > s sm2 ^ s + V -p> sV p> s cos2 Pp> s 
+V^ s U PiS e iSp < s cos p P:S sin [3 PjS (B.8) 
+U p , s V_ P)S e- i5p - a cos/3 Pifl sin/? P)S } . 

This corresponds to replacing 

V- P , s (t) => cos((3 p , s )V^ s (t) + e iS *»° sm(f3 p , s )U p , s (t) (B.9) 

in the original expression Q2.25| ). With the same substitution one obtains the 
expectation values for the energy and pressure. For the particle number the 
substitution is made in the Bogoliubov coefficients, Eqs. (|3.9|) and ( |3.15|) . 
As in the case of scalar fluctuations f3"7|] , we do not need the phase # P)S as 
long as we restrict the initial condition for the classical field to 0(0) = 0. 
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Figure Captions 

Fig. 1: <f)(t) for M = 1, g = 2, A = 1, and 0(0) = 0.6, 0(0) = 0. 

Fig. 2: Classical energy (dashed line), fluctuation energy (dotted line), 
and total energy (solid line) for the same parameters as in Fig. 1. 

Fig. 3: Total pressure for the same parameters as in Fig. 1. 

Fig. 4: Total particle number M\(t) (solid line) and Mo{t) (dotted line) 
for the same parameters as in Fig. 1. 

Fig. 5: Maximal occupation number J\f env as a function of E for the same 
parameters as in Fig. 1. 

Fig. 6: Exact quantum evolution (dashed line) and linearized evolution 
(solid line) of 0(t) for M = 1, g = 2, A = 1, and 0(0) = 0.01, 0(0) = 0. 

Fig. 7: Total pressure for the exact evolution for the same parameters 
as in Fig. 6. 

Fig. 8: Maximal occupation number J\f em as a function of E for the same 
parameters as in Fig. 6. 

Fig. 9: Total particle number Af±(t) for the same case as in Fig. 6. 

Fig. 10: Exact quantum evolution (dashed line) and linearized evolution 
(solid line) of 0(t) for M = 1, g = 1, A = and 0(0) = 0.1, 0(0) = 0. 

Fig. 11: 0(t) for M = 1, g = 2, A = 1 and 0(0) = 2, 0(0) = 

Fig. 12: Classical energy (dashed line), fluctuation energy (dotted line), 
and total energy (solid line) for the same parameters as in Fig. 11. 

Fig. 13: Maximal occupation number M eiw as a function of E for the 
same parameters as in Fig. 11. 

Fig. 14: Maximal occupation number A/" cnv for m(t) = mo cost with mo = 
0.1 (dashed line),m = 0.5 (solid line), and m = 1 (dotted line), as a func- 
tion of E . 
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Fig. 15: <f>(t) including back-reaction of both fermionic and scalar fluc- 
tuations for M = 1, g = 2, A = 4, and 0(0) = 1, 0(0) = 0. 

Fig. 16: (p(t) in the absence of scalar fluctuations for the same parame- 
ters as in Fig 15. 

Fig. 17: (j)(t) in the absence of fermion fluctuations, for the same parameters 
as in Fig. 15. 

Fig. 18: Fluctuation energies of fermions for evolution without scalar fluctu- 
ations (solid line); fluctuation energies of fermions (dashed line) and bosons 
(dotted line) for combined evolution. 
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